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Comment on "Fock-Darwin States of Dirac 
Electrons in Graphene-Based Artificial Atoms" 

The authors of Ref.pJ have computed eigenenergies 
and eigenstates of a graphene dot in the presence of 
two-dimensional (2d) harmonic potential and a magnetic 
field by considering only basis states with energies larger 
than or equal to zero. We show that the wavefunctions 
of the nearly chiral Landau level (LL) obtained using 
this approximate method are inaccurate, and violate a 
Hellman-Feynman theorem. Corrections lead to signifi- 
cant changes of optical strengths involving nearly chiral 
states. Moreover, additional basis states lead to new en- 
ergy levels in the energy spectrum. 

In Ref.[l| the following type of basis states of the 
Hamiltonian matrix are used with energies larger than 
or equal to zero 
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The energies of these states are E n = sgn(n)^^y / 2|n[ 
(vp is the Dirac velocity and I is the magnetic length). 
First order perturbation theory gives the perturbed wave- 
functions of nearly chiral LL states with n — 0: IV'om) ~ 
\ipo,m) + Ci|^i,m+i) + C_i|i/'-i,m+i>- These states have 
all the same angular momentum J = —m — 1/2. We 
find, using C_i = —Ci, that the expectation value of 
the azimuthal velocity (vg) is two times bigger when 
both positive and negative energy basis states are in- 
cluded than when only positive energy basis states are 


included in the perturbation series: (V'oml^slV'o 
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Our diagonalization of 


Cl(0o,m|e J \(f>0,m+l 

the Hamiltonian matrix confirms the simple perturbative 
result above. Figflja) shows that the difference in (vg) is 
about a factor of two. 

The azimuthal velocity (vg) of states in a quantum dot 
with V(r) = 4«T 2 should approach the velocity along 
a Hall bar[2j in the Id harmonic potential ^kx 2 . It 
is easy to show that the graphene edge electrons of a 
Hall bar obey the Hellman-Feynman theorem ^ = 
(ipn,k\ v y\' l l ; n,k}, where the y-component of wavevector is 
denoted by k and the energy of the n'th LL by e n (k). 
Figflja) displays the values of (vg) for nearly chiral states 
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of the 


of a dot, and we see that they converge to 
Hall bar (the relevant quantum numbers of the two sys- 
tems are related by kl — y/2(m + l)[2j). When negative 
energy basis states are excluded in the Hamiltonian ma- 
trix the value of the azimuthal velocity (vg) is reduced 
significantly while the value 4 d€ g^ does not change as 
much. The Hellman-Feynman theorem is thus violated 


in this case. 

Moreover, numerous new states appear in the energy 
spectrum, for example, between the energies of (0, 0) and 
(0, 1), as shown in FigQJb). Note also there are negative 
energy levels, which are absent in Ref.[l[. Some of these 
new levels give rise to additional optical transitions with 
significant strength. The appearance of additional lev- 
els with positive energies is due to negative energy ba- 
sis states that get perturbed by the parabolic potential 
with their energies shifted upward. In addition, we find 
that the strength of the optical transition (0,0) —> (0, 1) 
at B=2T with a = 0.925 and polarization along x-axis, 
shown in Fig. 3 of Chen et al.Q, is about 2.6 times larger 
than the correct value obtained by including both posi- 
tive and negative energy basis states. Physical origin of 
this difference is anticrossing of (0, 0) with other levels in 
the interval B < 4T, see Figfjjb). The states (1,0) and 
(0, 1) also anticross below 2T with new levels arising from 
the presence of negative energy basis states. This effect is 
missing in Ref.[l|, and leads to additional modifications 
of optical strength. 



FIG. 1: (a) Velocity (v) of the eigenstate with n = 0. Here 
a = = 0.2. Values of velocities using 21 ba- 

is for 2d harmonic potential (circles), and 


sis states: 

\ gl! 3fc fc ^ ijr are f° r Id and 2d harmonic potential systems (tri- 
angles), is for 2d harmonic potential computed exclud- 
ing negative energy basis states (squares). 11 basis states are 
included, (b) Energy levels in the Hilbert subspaces with an- 
gular momentum J = —1/2 (filled circles) and J = 1/2 (open 
circles) are shown together. The number basis states used is 
30. Allowed optical transitions obey AJ = ±1. 
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